Regge poles of the Schwarzschild black hole: a WKB approach 
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We provide simple and accurate analytical expressions for the Regge poles of the Schwarzschild 
black hole. This is achieved by using third-order WKB approximations to solve the radial wave 
equations for spins 0, 1, and 2. These results permit us to obtain analytically the dispersion relation 
and the damping of the "surface waves" lying close to the photon sphere of the Schwarzschild black 
hole and which generate the weakly damped quasinormal modes of its spectrum. Our results could 
be helpful in order to simplify considerably the description of wave scattering from the Schwarzschild 
black hole as well as the analysis of the gravitational radiation created in many black hole processes. 
Furthermore, the existence of nonlinear dispersion relations for the photons propagating close to 
the photon sphere could have also important consequences in the context of gravitational lensing. 

PACS numbers: 04.70.-s, 04.30.-w 



I. INTRODUCTION 



Since the sixties, mainly under the impetus of Nussen- 
zveig in electromagnetism and Regge in quantum physics, 
asymptotic/semiclassical techniques which use analytic 
continuation of partial-wave representations have been 
developed to understand scattering theory and the as- 
sociated resonance phenomenons. Together these tech- 
niques form the complex angular momentum (CAM) 
method. The CAM method originates from the pioneer- 
ing work of Watson [l[ dealing with the propagation and 
diffraction of radio waves around the earth (see also the 
modification proposed by Sommerfeld Q of the Watson 
approach). Today, the CAM method has been success- 
fully introduced in various domains of physics. For re- 
views of the CAM method we refer to the monographs 
of Newton Nussenzveig [1], Grandy [1], Collins @ 
and Gribov [7| as well as to references therein for various 
applications in quantum mechanics, nuclear physics, par- 
ticle and high-energy physics, electromagnetism, optics, 
and seismology. 

The CAM method is based on the dual structure of 
the S'-matrix associated with a given scattering problem. 
Indeed, the S'-matrix is a function of both the frequency 
w and the angular momentum index £. It can be analyt- 
ically extended into the complex uj plane as well as into 
the complex £ plane (CAM plane). Its poles lying in the 
fourth quadrant of the complex uj plane are the complex 
frequencies of the resonant modes. In other words, the 
behavior of the ^-matrix in the complex uj plane permits 
us to investigate resonance phenomenons. The struc- 
ture of the iS-matrix in the complex £ plane allows us 
- by using integration contour deformations, Cauchy's 
Theorem, and asymptotic analysis ~ to provide a semi- 
classical description of scattering in terms of geometric 



and diffractive rays. In that context, the poles of the S'- 
matrix lying in the first quadrant of the CAM plane (the 
so-called Regge poles) are associated with diffraction. Of 
course, when a connection between these two faces of the 
S'-matrix can be established, resonance aspects of scat- 
tering are then semiclassically interpreted. 

More precisely, for spherically symmetric scattering 
problems, the connection between the resonant approach 
of scattering and its Regge pole interpretation is a con- 
sequence of the following considerations: 

i) For a given value £ G N of the angular momentum 
index, the resonances, i.e., the complex poles of the ma- 



trix element Si{uj), are of the form uj£, 



- ir£„/2, 



where n G N and, in the immediate neighborhood of uj^p, 
Si{uj) has the Brcit-Wigncr form, i.e.. 



Siiuj) cx 



(1) 



ii) For a given value cj > of the frequency, the 
Regge poles, i.e., the complex poles of the matrix el- 
ement Sa_i/2(w), can be expressed as function of the 
form A„(a;). The curves traced out in the CAM plane 
by the Regge poles as a function of the frequency lu are 
the so-called Regge trajectories. They permit us to inter- 
pret Regge poles in terms of diffractive rays or "surface 
waves": Re A„(cj) provides the dispersion relation for the 
nth "surface wave", while ImA„(w) corresponds to its 
damping. 

(iii) From the Regge trajectories, we can semiclas- 
sically construct the resonance spectrum. There exists 
a first semiclassical formula (a Bohr- Sommerfeld- type 
quantization condition) which provides the location of 

the excitation frequencies uj^^ of the resonances gener- 
ated by nth "surface wave" : 



RcK{ujfj) =£+1/2, £ 



e N. 



(2) 
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There exists a second semiclassical formula which pro- 
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vides the widths of these resonances 



Im Xn{uj)[d/ du Re A„(a;)] 



~ ^ '[d/da;ReA„(a;)]2 



[d/dw Im A„(aj)] 



(3) 



and which reduces, in the frequency range where the con- 
dition jd/dwRe A„(cj)| ^ |(i/(iajImA„(w)| is satisfied, to 



r^n _ ImA„(^) 
2 d/dw ReA„(cij) 



(4) 



Recently, by using CAM techniques or in other words 
the Regge pole machinery, we have established that the 
weakly damped quasinormal mode (QNM) complex fre- 
quencies of the Schwarzschild black hole of mass M are 
Breit-Wigner-type resonances generated by a family of 
"surface waves" lying on its photon sphere at r = 3M 
3. More precisely, by noting that each "surface wave" 
is associated with a Reggepole of the ^-matrix of the 
Schwarzschild black hole |8l-[lO|. we have been able to 
numerically construct the spectrum of the QNM com- 
plex frequencies from the associated Regge trajectories. 
In fact, in this way, from CAM methods, we have given 
meaning to an appealing and physically intuitive inter- 
pretation of the Schwarzschild black hole QNMs sug- 
gested by Goebel in 1972 i.e., that they could be 



interpreted in terms of gravitational waves in spiral or- 
bits close to the unstable circular photon orbit at r = 3M 
which decay by radiating away energy. It should be noted 
that alternative implementations of the Goebel interpre- 
tation have been developed (see Refs. [l^-[I3) but they 
are mainly based on the concepts of geodesic and bundle 
of geometrical rays while our CAM analysis is based on 
wave physics. 

In Sec. 2 of this short note, we provide analytical ex- 
pressions for the Regge poles corresponding to the scalar 
field (spin 0), the electromagnetic field (spin 1) and the 
gravitational perturbations (spin 2) propagating on the 
Schwarzschild background. This is achieved by using the 
WKB approach developed in the context of the determi- 
nation of the QNMs by Schutz and Will [II and by Will 
and Iyer [l^ [23| (see also Ref. [2l| for general aspects of 
WKB theory and for particular aspects connected with 
eigenvalue problems of the type considered here). Our 
WKB results permit us to describe very accurately the 
Regge trajectories of the Schwarzschild black hole and 
to recover from the semiclassical formulas ^ and ^ 
well-known analytical expressions for the QNM complex 
frequencies. It should be noted that in Ref. @, from 
numerical considerations, we have provided a coarse for- 
mula for the Regge poles of the Schwarzschild black hole. 
It has been sufhcient in order to obtain a "surface wave" 
interpretation of the QNM complex frequencies but it 
does not take into account the spin dependent aspects of 
the problem as well as the fact that (i) the dispersion re- 
lation of the nth "surface wave" is nonlinear and depends 
on the index n and that (ii) the damping of the nth "sur- 
face wave" is frequency dependent. In the present pa- 
per, we have derived analytical expressions for the Regge 



poles of the Schwarzschild black hole from mathemati- 
cal considerations. The results displayed in Ref. Q are 
the leading-order counterpart for very high frequencies 
of our new results which are much more accurate and 
which now permit us to describe precisely the behavior 
of the Regge poles and of the associated "surface waves" 
in a very large range of frequencies. In Sec. 3, we briefly 
conclude by considering some possible consequences and 
applications of our results. 



II. REGGE POLES, REGGE TRAJECTORIES, 
AND WKB RESULTS 



The wave equations for the scalar field, for the elec- 
tromagnetic field, and for the gravitational perturbations 
propagating on the Schwarzschild black hole reduce, after 
separation of variables, to the Regge- Wheeler equation (a 
one-dimensional Schrodinger equation) 



dr'i 



where the effective potential V^{r) is given by 



Vt{r) = 



2M 



£{£+!) 2(l-s2)Af' 



(5) 



(6) 



Here, s denotes the spin of the field considered, £ e N 
is the ordinary angular momentum index, while r and 
r,, — r + 2Mln {r/2M — 1) -|- const are, respectively, the 
standard radial Schwarzschild coordinate and the Regge- 
Wheeler tortoise coordinate, and we have furthermore 
assumed a harmonic time dependence e~*'^*. 

For a given angular momentum index £, the ^-matrix 
element Si{uj) is obtained by seeking the solutions of the 
Regge- Wheeler equation ([5]) which have a purely ingoing 
behavior at the event horizon r ~ 2M, i.e. which satisfy 



$^(r) - 

— f — CO 

and which, at spatial infinity r 
totic behavior of the form 



(7) 



-00, have an asymp- 



$^(r) 



1 



T,(w)' 



(8) 

For w > 0, the Regge poles of the ^-matrix are 
the complex values A for which, both S'a-i/2('^) a-nd 
T\^i/2{^) have a simple pole but S'a- 1/2 ('^)/r'A- 1/2(1^) 
is regular [8]. They constitute a family A„(aj) with the in- 
dex n = 1,2,... permitting us to distinguish between the 
different poles. The corresponding modes are the solu- 
tions of the radial wave equation ([5]) which are purely out- 
going at infinity and purely ingoing at the horizon. They 
are the "Regge modes" of the Schwarzschild black hole. 
It should be noted that the QNMs of the Schwarzschild 
black hole and the associated complex frequencies are de- 
fined in an analogous way with moreover a fundamental 
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FIG. 1: Regge pole distributions for uj = 1 and uj = 2 {2M — 
1). Scalar and electromagnetic fields as well as gravitational 
perturbations are considered. 
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FIG. 2: Regge trajectories for the scalar field. The Regge 
poles A„(cj) (n = 1, 2, 3) are followed for cj = ^> 5 (2M = 1). 



difference. Indeed, they are obtained for fixed values of 
£ e N and exist only for complex values of uj. 

In Ref . ^ , we have noted that the Regge poles can be 
determined by using, mutatis mutandis, the method de- 
veloped by Leaver [22| for the search of the QNM complex 
frequencies and we have furthermore numerically imple- 
mented this method from the Hill determinant approach 
of Majumdar and Panchapakesan [2^. For the Regge 
poles lying close to the real axis of the CAM plane, the 
"exact" results we can obtain numerically are accurate 
ones. Figure [T] exhibits the distribution of Regge poles 
for two different values of the frequency a;. On this fig- 
ure, the splitting of the Regge poles with the spin s and 
the index n as well as the migration of the first Regge 
poles corresponding to spins 0, 1 and 2 can be observed. 
This migration can be more precisely described by study- 
ing the Regge trajectories, i.e., the curves traced out in 
the CAM plane by the functions A„(a;) for uj S [0, -|-oo[. 
In Figs, imi] we have displayed the Regge trajectories for 
the first three Regge poles associated with the scalar field, 
the electromagnetic field and the gravitational perturba- 
tions. As we have already noted in Ref. [8], Re A„(w) and 
Im A„ (uj) respectively denote the azimuthal propagation 



constant (i.e., the dispersion relation) and the damping 
constant of the nth "surface wave" lying on the photon 
sphere. So Figs. permit us to note that, at first sight, 
(i) the global behavior of these "surface waves" is rather 
independent of the spin considered (such result is only 
valid for high frequencies) and (ii) for a given spin, the 
dispersion relation seems to be independent of the "sur- 
face wave" considered (such result is also only valid for 
high frequencies) and the "surface wave" corresponding 
to the first Regge pole has the weakest attenuation and 
therefore must play the most important role in scatter- 
ing and in radiative processes and must contribute sig- 
nificantly to the resonance mechanism. 

The WKB approach developed in the context of the 
determination of the QNMs by Schutz and Will [fSi] and 
by Will and Iyer [H,!!^ (see also Ref. ^l|) can be easily 
adapted in order to obtain analytical approximations for 
the Regge poles. Indeed, as we have already noted, Regge 
modes and QNMs both satisfy the same wave equation 
and the same boundary conditions but with different in- 
terpretations for the angular momentum index and the 
frequency. By using third-order WKB approximations 
[l^, [1^ for the Regge modes, we can find that the Regge 
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Regge trajectories for gravitational perturbations (2M=1) 
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FIG. 3: Regge trajectories for the electromagnetic field. The 
Regge poles A„(aj) (n = 1,2,3) are followed for lj = — s- 5 
(2M = 1). 



FIG. 4: Regge trajectories for gravitational perturbations. 
The Regge poles \n{oj) (n = 1, 2, 3) are followed for a; = — 5- 
5 (2M = 1). 



and 



poles A are the complex solutions of the equation 



V^{\) + [~2V^^>{\)] A(A,n) 
ia{n) [^2V^''\X)]'^[l+Ti{\,n)] (9) 



w^ith > and n — 1,2,3,.... Here 



A(A,n) = 



1/2 



288 [v^-\x)^ 



(7 + 60 a{nf 



(10a) 



-2V^'\X)] 
6912 [v^^)^x)) 



(77+ 188a(n)2) 



1 
384 



.(4) 



+ 



1 





'Vo^'\x)V,^'\xy 



288 V [V,^'\X)] 

1 fv,^'\x)' 



{19 + 28 a{nf) 



(5 + 4a(7i)2 



(10b) 



288 \Vo^'\x), 
In Eqs. (0) and ([TU]), we have introduced the notations 

Q;(n)=n-l/2 (11a) 



5 



and, for p G N, 



dp 



dr.. 



(lib) 



with (r*)Q which denotes the maximum of the function 
Vx_i/2{r*)- We can solve Eq. © by assuming |A| ^ 1 
as weU as Re A ^ ImA. Its solutions constitute a family 
A„(a;) with n = 1, 2, 3, . . . , and we have the approxima- 
tion 



A„(c^) « [27{Muj)^ + an + 2(n - l/2)2e„(u;)]'/' 



+^ (n- l/2)[l + e„(tj)] 



where 



and 



with 



12 



(n- 1/2)2 + 



115 

144 



1 + 



27(Afa;)2 + a„ + (n- 1/2)2 



, 4 

o„ = — 
9 



305 
1728 



(n- 1/2)2 + 



5555 
6912 



- 1 + .5^ 



(12) 



(13) 



(14) 



(15) 



Equation (jl2l) is the main result of our paper. It also 
provides expressions for the dispersion relation and the 
damping of the "surface waves" lying on/close the pho- 
tosphere of the Schwarzschild black hole. It is important 
to note that the perturbation parameter of the WKB 
method leading to Eq. ^ is not 1/uj with w — +oo but, 
in some sense, the distance e-^^^ (with formally e — > 0) 
between the turning points of Eq. ^ [i.e., the roots of 
Q{r*) = uP' — V^(r*)] and the location of the peak of 
Q(r*) (see also Refs. |18l-t20|). As a consequence, Eq. (fT2|) 
is valid in a very large range of frequency values and not 
only for u) — > +oo. Of course, in order to solve Eq. © 
and to obtain the expression (|12p. we have furthermore 
assumed that |A| ^ 1 as well as Re A ^ ImA. As a 
consequence, we can expect the accuracy to decrease for 
very low frequencies and for higher Regge poles. 

Tables Hl lIIII present a sample of Regge pole positions 
calculated from the WKB resuh (HH). The WKB resuhs 
are better for high frequencies and for the first Regge 
pole. However, a comparison between the "exact" and 
the WKB results shows a good agreement even for low 
frequencies and for the second and third Regge poles. 
Formula (jl2[) permit us also to describe very accurately 
the Regge trajectories of the Schwarzschild black hole 
[see Fig. [S] where, for clarity, we have considered the 
Regge trajectories in the form ImA„(aj) as a function 
of ReA„(w)]. It should be also noted that, in Ref. 0], 
Glampedakis and Andersson have developed a "quick and 
dirty" method which provides very accurate results for 
the Regge poles (see Table 4 in their article where the 
first three Regge poles associated with the scalar field 



TABLE I: A sample of Regge pole values for the scalar field 
(2M = 1): Exact versus WKB results. 
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are given). Their numerical method is certainly better 
than the WKB method we have considered. But it is 
only a purely numerical method used in order to approx- 
imate the exact results. Our WKB approach permits us 
to provide a simple and accurate analytical expression 
for the Regge poles. It would be also possible to improve 
considerably our WKB results by adapting the higher or- 
der WKB approach developed by Konoplya in Ref. '255 ■ 
But we believe that the Konoplya approach could not 
provide a formula as simple as (|12p. 

It is moreover possible to simplify (|12p and to write 



A„(a;) = 



+i {n - 1/2) 



ISMw 
b„ 



27(Ma;)2 



O 



(16) 



which is much more elegant. Of course, this alterna- 
tive formula provides accurate results only at large lu 
and must not be used at low frequencies (see Fig. |6]). It 
provides very simple expressions for the dispersion rela- 
tion and the damping of the "surface waves" lying on 
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TABLE II: A sample of Regge pole values for the electromag- 
netic field (2Af = 1): Exact versus WKB results. 
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the photosphere of the Schwarzschild black hole. The 
leading-order terms of Eq. ([T6| correspond to the ex- 
pression (19) of Ref. [1]. Of course, this expression is too 
coarse to take into account precisely the spin and index 
dependence which appear clearly in Fig. [1] and, further- 
more, it does not emphasize the nonlinear behavior of the 
dispersion relations as well as the nonconstant behavior 
of the damping (see Figs. [Hill). As a consequence, the 
numerical calculations which use it are not accurate for 
low frequencies uj, for high indices n and for high spins s. 
For example, for 2M = 1, w = 1, n = 3, and s = 0, the 
formula (19) of Ref. [1] leads to important errors about 
8.6 % for the real part of the Regge pole and —4.5 % for 
its imaginary part. For s = 1 the errors are, respectively, 
about 10.0 % and —7.3 %, and for s = 2 they are about 
14.8 % and -16.2 %. 

In a recent paper [26| written after the first version 
of the present one, Dolan and Ottcwill have greatly im- 
proved our high frequency approximation (|16p by pro- 
viding higher orders in 1/uj. But, as they have noted, 
their expansion is not as accurate as our formula (jl2p in 
the low-w regime. This reinforces again our view that 
Eq. (|12p is a very interesting result due to its large range 



TABLE III: A sample of Regge pole values for gravitational 
perturbations (2M = 1): Exact versus WKB results. 
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of validity. 

In Ref. [1] , we have proven that the "surface waves" as- 
sociated with the family of Regge poles propagate on the 
photon sphere of the Schwarzschild black hole at r = 3M. 
We have established this connection by using the leading- 
order behavior of (fTB)) . It is therefore a semiclassical in- 
terpretation which is formally valid for w — >■ +oo. In this 
context, we recall (see Ref. Q) that the resonant part of 
the form factor of the Schwarzschild black hole is a su- 
perposition of terms like exp[i{Xn{uj)0 — wf)] which can 
be also written in the form 

exp[-/c;'(c^)/:] exp[z(fc;(c^)/: - ujt)\ (17) 

with £ = (3Af)6' which denotes the arc length taken on 
the photon sphere. Here k'^^iyi) = ReA„((u;)/(3A/) rep- 
resents the wavenumber of the nth "surface wave" lying 
on the photon sphere, or in other terms its dispersion 
relation, while fc"(i^) = Im A„(cj)/(3M) is its damping 
constant. To the leading order, the dispersion relation is 
linear but if we go beyond this order, it is clearly nonlin- 
ear [see Eqs. and (|16p ]. and this last result could have 
important consequences in strong gravitational lensing. 
But the previous analysis could be also inexact: indeed. 
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Regge trajectories for gravitational perturbations (2M=1) 
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FIG. 5: Regge trajectories spins 0, 1 and 2: Exact versus 
WKB calculations. For the Regge poles An(i^) (n = 1,2,3) 
followed for a; = — >■ 8 (2A/ = 1), ImA„(cj) is expressed as a 
function of Re A„(cj). 



if we need to take into account the finite value of the pho- 
ton frequency, the nonlinear corrections could also induce 
a shift for the location of the unstable circular orbit of 
the photon and it would be then necessary to slightly 
modify (jf 7p and more particularly the expression of the 
arc length £. 

Finally, let us note that our WKB results permit us to 
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FIG. 6: Regge trajectories (spin 2, n = 1): Exact versus 
WKB calculations. The Regge pole Ai(aj) followed for uj = 
0.f2 1.83 (2M = 1). 



recover from the semiclassical formulas ([2]) and (jlj well- 
known analytical expressions for the QNM complex fre- 
quencies. Indeed, by inserting (IT^ into Eqs. ^ and (jH), 
we obtain the well-known large I behaviors [see Eq. (3.1) 
ofRef. ^0^] 
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(19) 

Here, it is worth recalling that the leading order of (|18ap 
can be expressed in terms of the angular velocity at the 
unstable circular null geodesic, i.e. at r = 3M, while 
the leading order of (|18b|) is related to the corresponding 
Lyapunov exponent (see Ref. 17]). 



III. CONCLUSION 

The success of the CAM method in physics is due to its 
ability to provide a clear description of a given scattering 
problem by extracting the physical information (linked 
to the geometrical and difFractive aspects of the scatter- 
ing process) which is hidden in partial-wave representa- 
tions. In the context of black hole physics, the CAM 
method has permitted us to clearly interpret the exis- 
tence of the weakly damped QNM complex frequencies of 
the Schwarzschild black hole in terms of "surface waves" 
propagating on/close to its photon sphere [8l-[l0| and, in 
the present paper, we have derived analytical expressions 



8 



for the dispersion relation and the damping of these "sur- 
face waves" [see Eqs. and p^ ] or, equivalently, for 
the location in the CAM plane of the Regge poles of the 
Schwarzschild black hole associated with these "surface 
waves" . 

The analytical expressions we have obtained for the 
Regge poles of the Schwarzschild black hole may be help- 
ful in order to simplify considerably (i) the description 
of wave scattering from this gravitational background as 
well as (ii) the description of the gravitational radiation 
created in many black hole processes involving this gravi- 
tational background. In the former context, the approach 
initiated by Andersson and Thylwe d, [13] is certainly the 
most promising way. In the latter one, a formalism based 
on the Regge pole machinery has not yet be developed 
and all the work remains to be done. 

Our results could have also very important conse- 
quences and applications in the context of gravitational 
lensing. Indeed, until now the gravitational lens effects 
due to the Schwarzschild black hole have been only con- 
sidered in the framework of geometrical optics where light 
propagation is described in terms of null geodesies. In 
this framework, photon dispersion is completely irrele- 
vant. But, in the present paper, we have established 
from wave theoretical considerations the existence of dis- 
persion relations for the photons propagating close to 
the photon sphere at r = 3M and it is well-known that 
this "surface plays" a crucial role in the context of the 
strong-deflection gravitational lensing (see, for example, 
Ref. 27] and references therein). As a consequence, it 
seems to us very important to reexamine the theoreti- 



cal analysis of gravitational lensing in the light of our 
new results. But even without doing such a work, we 
can already envision some of their possible observational 
consequences such as iridescence phenomenons for the 
images produced by the Schwarzschild lens, time delay 
for photons of different frequencies as well as a hyperfine 
structure in the system of images (due to the existence 
of an infinity of Regge poles) . 

In the sixties, high-energy physics was mainly based 
on the Regge pole approach of scattering processes (see, 
for example, Refs. [y, 0] and references therein). Of 
course, after the discovery of quarks and the develop- 
ment of Yang-Mills theories, it has taken completely dif- 
ferent paths. So, it is amazing to encounter again the 
Regge pole concept but now in the context of gravita- 
tion physics (see also Ref. [11] where Regge poles seem 
to play an important role in the context of the AdS/CFT 
correspondence). It is then rather natural to wonder if 
there could be something very fundamental behind this 
which, in particular, could be helpful in order to under- 
stand black holes and gravitation at the quantum level. 
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